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Abstract 

Following the paradigm of Boltzmann-BBGKY we propose a correlation entropy (of the nth 
order) for an interacting quantum field, obtained by 'slaving' (truncation with causal factor- 
ization) of the higher (n+1 th) order correlation functions in the Schwinger-Dyson system of 
equations. This renders an otherwise closed system effectively open where dissipation arises. 
The concept of correlation entropy is useful for addressing issues related to thermalization. As 
a small yet important step in that direction we prove an H-theorem for the correlation entropy 
of a quantum mechanical O(N) model with a Closed Time Path Two Particle Irreducible Effec- 
tive Action at the level of Next-to-Leading-Order large N approximation. This model may be 
regarded as a field theory in space dimensions. 

1 Introduction 

The goal of this paper is two-fold: the proposal of a correlation entropy for an interacting quantum 
field, and the proof of a H-theorem for the quantum mechanical O(N) model ^ |2j, which may 
be regarded as a field theory in zero space dimensions. For the former, we follow the paradigm 
of Boltzmann-BBGKY [3| and propose a correlation entropy (of the nth order) for an interacting 
quantum field [3J 03 |S| obtained by 'slaving' the higher (n+1 th) order correlation functions in 
the Schwinger-Dyson system of equations (See [7J |S] and below). We then derive the closed time 
path (CTP) [S] two particle irreducible (2PI) ^U] effective action (EA) for this model up to the 
next-to-leading order (NLO) in 1/N As a useful step towards addressing issues related to the 
thermalization of interacting quantum fields, we prove an H-theorem for the correlation entropy 
of a quantum mechanical O(N) model at the CTP 2PI NLO level. We give an introduction to 
the correlation hierarchy and the thermalization issue in this section and a short description of 
different definitions for the entropy of interacting quantum fields in the next section, followed by a 
derivation of the H-theorem given in the ensuing sections. 1 

*Email: calzetta@df.uba.ar 
^Email: hub@physics.umd.edu 

1 The following descriptive summary in this and the next section is adapted from Chapter 5 of |12|. where a 
preliminary attempt to construct a correlation entropy was made but the proof of the H-theorem for an interacting 
quantum field did not materialize, and from |13|. 
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1.1 Boltzmann-BBGKY Hierarchy and Schwinger-Dyson Equations 

We begin with nonequilibrium statistical mechanics. As is well known, truncation of the BBGKY 
hierarchy at a finite nth order yields a closed system 2 . The system of equations describing the 
n-particle distribution functions are time reversal invariant. When a causal factorization condition 
such as Boltzmann's molecular chaos assumption is imposed - that is, assuming that the n+1 th 
correlation function can be factorized into a product of n th correlation functions initially but not 
finally (after the collisions) - time-irreversibility appears and an //-theorem obtains. This type 
of coarse graining of the hierarchy i.e., truncation plus causal factorization into a set of coupled 
equations for the n-particle distribution functions is called 'slaving' in the language of IB]- Slaving 
of the n+1 -particle distribution function renders an otherwise closed system (of the n+1 th order 
correlation functions) an effectively open system (of nth order correlation functions) and ushers in 
the appearance of dissipative dynamics UJ. Noise and fluctuations ^JE3 should also appear, as 
required by the fluctuation-dissipation relation (FDR) now manifesting not just for an open 
system near equilibrium (as depicted in linear response theory), but for an effectively open system 

A familiar example is the Boltzmann equation: it is not enough to truncate the BBGKY 
hierarchy to include only the one-particle and two-particle distribution functions. One needs to 
make an additional assumption that the two-particle distribution function factorizes into a product 
of single-particle distribution functions before the collisions, i.e., invoking the molecular chaos 
hypothesis. Of course the colliding partners become correlated after interaction. This causal 
distinction is the origin of the macroscopic arrow of time and the appearance of dissipation in the 
Boltzmann paradigm. FDR calls for a rightful place for fluctuations in the Boltzmann equation, as 
demonstrated in the derivation of a stochastic Boltzmann equation |171 118j . 

The situation in quantum kinetic field theory is completely analogous. By quantum kinetic field 
theory (see, e.g., [U 03 E] ) > we are referring to the hierarchy of coupled equations for the relativis- 
tic Wigner function and its higher-correlation analogs, obtainable from the variation of the master 
effective action [315] whose variation yields the Schwinger-Dyson equations. This is a quantum ana- 
logue of the BBGKY hierarchy, expressed in a representation convenient for distinguishing between 
microscopic (quantum field-theoretic) and macroscopic (transport and relaxation) phenomena 3 . 

One may choose to work with a truncation of the hierarchy of the Wigner function and its 
higher correlation analogs, or one may instead perform a slaving of, for example, the Wigner- 
transformed four-point function, which leads (within the context of perturbation theory) directly 
to the relativistic Boltzmann equation and the usual //-theorem. The truncation and subsequent 
slaving of the hierarchy within quantum kinetic field theory can be carried out at any desired order, 
as dictated by the initial conditions and relevant interactions. 

1.2 Thermalization in Interacting Quantum Fields 

The problem of thermalization in relativistic quantum fields has drawn much attention over time, 
both in the attempt to understand the origin of macroscopic irreversible behavior from microscopic 
theories, and for practical considerations of nonequilibrium quantum field processes in the early 

2 The truncation of the correlation hierarchy is not just an arbitrary mathematical procedure, it reflects the fact 
that in realistic conditions, measurement are of finite accuracy associated with limited resolution of the instruments. 
Thus the relevant physical degrees of freedom are often limited to the lower end of the hierarchy, namely, the mean 
field and two-point function. 

3 It should be pointed out that in order to identify the relativistic Wigner function with a distribution function for 
quasiparticles, one must show that the density matrix has decohered, and this is neither guaranteed nor required by 
the existence of a separation of macroscopic and microscopic time scales |19| . 



2 



Universe and in relativistic heavy ion collisions |2Uj . 

Let us begin with a discussion of the exact meaning of thermalization. In the strictest sense 
an isolated system depicted by quantum field theory undergoes unitary evolution and does not 
thermalize. However, one can still ask meaningful questions such as whether certain correlation 
functions may converge to their thermal forms in some well defined physical limit. We may call 
this a weak thermalization condition. 

In practical terms, even asking this type of questions may seem too ambitious. As described 
above, the system of equations for the correlation functions form an infinite hierarchy. When cer- 
tain causal factorizable boundary condition is mposed, the truncated subsystem will show signs of 
irreversibility and a tendency to equilibrate. For interacting fields we have to deal with a truncated 
Schwinger - Dyson hierarchy to make the analysis possible. Moreover, nontrivial (point) field theo- 
ries are plagued by divergences which can only be controlled by regularization and renormalization 
within some perturbative scheme. Therefore the question of even weak thermalization can make 
sense only within a chosen level of approximation, such as the nPI of the correlation hierarchy, the 
loop expansion, the 1/N expansion, or expansion in powers of a coupling constant, etc. 

There have been recent claims based on numerical evidence [211 122j that at the next-to-leading 
order (NLO) Large N approximation an interacting quantum field may show signs of thermalization. 
In a recent report ^H] we consider a 0(N) invariant scalar field of unbroken symmetry, develop the 
CTP 2PI EA 6 in powers of 1/N, retaining up to the next-to-leading order [O(l)] terms, and show 
that the only time translation invariant solutions are thermal. Our analytic result provides support 
for similar claims. Here, we present an alternative approach by defining a correlation entropy and 
showing the existence of an H-theorem in the simpler quantum mechanical O(N) model. 

2 Entropy of Interacting Quantum Fields 

Because entropy is of such fundamental importance, it is useful to discuss how entropy can be 
defined for the various ways of approximating the features and dynamics of a quantum field. First, 
for a unitarily evolving quantum field theory whose dynamics is a closed system (as opposed to an 
'effectively open system') and governed by the quantum Liouville equation, it is well known that 
the von Neumann entropy of the density matrix, 

S V N = -Tr[p(t)lnp(t)}, (1) 

is exactly conserved. If there is a justifiable separation of macroscopic and microscopic time scales, 
one can adopt the theoretical framework of quantum kinetic field theory. If one makes the as- 
sumption of factorization (equivalently, slaving of the Wigner-transformed four-point function), 
one obtains the relativistic Boltzmann equation in the binary collision approximation. The Boltz- 
mann entropy Sb defined in terms of the phase space distribution f(k, X) for quasiparticles can 
in this case be shown to satisfy a relativistic //-theorem |23| 16] . We want to generalize this to a 
correlation entropy for interacting quantum fields. 

However, in the case where there does not exist such a separation of time scales, how does 
one define the entropy of a quantum field? For nonperturbative truncations of the dynamics of 
interacting quantum fields, this is a nontrivial question |24j . Intuitively, one expects that any 
coarse graining which leads to an effectively open system with irreversible dynamics will also lead 
to the growth of entropy. Formally these operations can be systematically expressed in terms of the 
projection operator techniques [22]. A projection operator P projects out the irrelevant degrees 
of freedom (thus going over to an open system) from the total system described by the density 
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operator p, yielding the reduced density matrix p R 



p R (t) = Pp(t) 



(2) 



There exists a well-developed formalism for deriving the equation of motion of the relevant degrees 
of freedom, and in terms of it, the behavior of the coarse-grained entropy |251 13] 



which will in general not be conserved. The projection operator formalism can be used to express 
the slaving of higher correlation functions in the correlation hierarchy. From it one can define an 
entropy in effectively open systems (see, e.g., |26p. (So far it has only been implemented within 
the framework of perturbation theory.) Another equally powerful method adept to field theory is 
the Feynman- Vernon influence functional formalism |27l I15j which has been used to treat entropy 
in quantum open systems (see, e.g., 28 J. 

2.1 Entropy special to Choice of Basis or Representation 

We now consider the entropy functions for quantum fields, beginning with the simpler yet more 
subtle case of a free field. Historically this issue was, to our knowledge, first raised in the context 
of entropy generation from particle creation for a free quantum field in an expanding Universe |22] 
due to the parametric amplification of vacuum fluctuations. The focal point is a wave equation 
with a time-dependent natural frequency for the amplitude function of a normal mode. (The same 
condition arises for an interacting field (such as the A<3> 4 theory) in the Hartree-Fock approximation 
or the O(N) field theory at leading order in the large- N expansion jllj. 1 ) Since the underlying 
dynamics is clearly unitary and time-reversal invariant in this case, a suitable coarse graining 
leading to entropy growth is not trivially evident. Hu and Pavon [3U] first made the observation 
that a coarse graining is implicitly incorporated when one chooses to depict particle numbers in 
the n-particle Fock (or "N") representation or to depict phase coherence in the phase (or "P") 
representation. Various proposals for coarse graining the dynamics of parametric oscillators have 
followed EH ESI OS EH] • The language of squeezed states is particularly useful for describing 
entropy growth due to parametric particle creation [SHI EZ1 ESI HE] • For our purposes, the essential 
features of entropy growth due to parametric particle creation which distinguish it from correlational 
entropy growth (to be discussed below), are that parametric particle creation depends sensitively 
on the choice of representation for the state space of the parametric oscillators, and the specificity 
of the initial conditions. 

2.2 Entropy from Projecting out Irrelevant Variables 

In contrast to entropy growth resulting from parametric particle creation from the vacuum, entropy 
growth due to particle interactions in quantum field theory [23] has a very different physical origin. 
A coarse graining scheme was proposed by Hu and Kandrup [21] for these processes. Expressing an 
interacting quantum field in terms of a collection of coupled parametric oscillators, their proposal 
is to define a reduced density matrix by projecting the full density operator onto each oscillator's 
single-oscillator Hilbert space in turn, 



and defining the reduced density operator as the tensor product n of the projected single-oscillator 



Scg 



Tr[p R (t) In p R (t)}, 



(3) 



7 (fc) = Trp^p, 



(4) 



density operators 7(/c), 



PR = %7(fe). 



(5) 
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The coarse-grained (Hu-Kandrup) entropy is then just given by Eq. (J2J), from which we obtain 



It is interesting to observe that for a spatially translation-invariant density matrix for a quantum 
field theory which is Gaussian in the position basis, this entropy is just the von Neumann entropy 
of the full density matrix, because the spatially translation-invariant Gaussian density matrix sep- 
arates into a product over density submatrices for each k oscillator. This projection (Hu-Kandrup) 
coarse graining, like the correlation-hierarchy (Calzetta-Hu) coarse graining scheme, does not choose 
or depend on a particular representation for the single oscillator Hilbert space. It is sensitive to 
the establishment of correlations through the explicit couplings. 

2.3 Entropy from Slaving the Higher Correlations: The CTP 2PI effective ao 



A general procedure has been presented for obtaining coupled equations for correlation functions at 
any order / in the correlation hierarchy, which involves a truncation of the master effective action 
at a finite order in the loop expansion |Hj . By working with an I loop-order truncation of the 
master effective action, one obtains a closed, time-reversal invariant set of coupled equations for 
the first I + 1 correlation functions, <f>, G, C3, . . . , C/+i. In general, the equation of motion for the 
highest order correlation function will be linear, and thus can be formally solved using Green's 
function methods. The existence of a unique solution depends on supplying causal boundary 
conditions. When the resulting solution for the highest correlation function is back-substituted 
into the evolution equations for the other lower-order correlation functions, the resulting dynamics 
is not time-reversal invariant, and generically dissipative. Thus, as was described before, with 
the slaving of the higher-order (Wigner-transformed) correlation function in quantum kinetic field 
theory, we have rendered a closed system (the truncated equations for correlation functions) into 
an effectively open system. In addition to dissipation, one expects that an effectively open system 
will manifest noise/fluctuations, as shown in |H] for the case of the slaving of the four-point 
function to the two-point function in the symmetry-unbroken A$ 4 field theory. Thus a framework 
exists for exploring irreversibility and fluctuations within the context of an unitary quantum field 
theory, using the truncation and slaving of the correlation hierarchy. The effectively open system 
framework is useful for precisely those situations, where a separation of macroscopic and microscopic 
time scales (which would permit an effective kinetic theory description) does not exist, such as is 
encountered in the thermalization issue. 

While it is certainly not the only coarse graining scheme which could be applied to an interact- 
ing quantum field, the slaving of higher correlation functions to lower-order correlation functions 
within a particular truncation of the correlation hierarchy, as a particular coarse graining method, 
has several important benefits. It can be implemented in a truly nonperturbative fashion. This 
necessitates a nonperturbative resummation of daisy graphs, which can be incorporated in the 
truncation/slaving of the correlation hierarchy in a natural way. 

In this study, we are interested in the growth of entropy due to the coarse graining of the 
correlation hierarchy by slaving a higher correlation function. The simplest nonperturbative trun- 
cation of the Schwinger-Dyson equations for the A<J> 4 field theory which contains the time-dependent 
Hartree-Fock approximation is the two-loop truncation of the master effective action, in which only 
the mean field cj>, the two-point function G, the three-point function C3 are dynamical. All higher 
order correlation functions obey algebraic constraints, and can thus be expressed in terms of the 
three dynamical correlation functions. 




(6) 
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While this truncation of the Schwinger-Dyson equations is well-defined and could in principle be 
solved, it is disadvantageous because, as stated above, without some coarse graining, the system will 
not manifest irreversibility and will not equilibrate. Therefore we slave the three-point function to 
the mean field and two-point function, and thus arrive at an effectively open system. In principle, a 
systematic analysis of the coarse-grained dynamics of the mean field and two-point function should 
include stochasticity |71|S], but noise is not our primary concern here. 

3 The H theorem for the quantum mechanical O (N) model 
3.1 The Next-to-Leading order Large N Approximation 

The number iV of replicas of essentially identical fields (like the N scalar fields in an O(N) invariant 
theory, or the N 2 — l gauge fields in a SU(N) invariant non-abelian gauge theory) suggests using 1/JV 
as a natural small parameter, with a well defined physical meaning. Unlike coupling constants, this 
is not subjected to renormalization or radiative corrections. By ordering the perturbative expansion 
in powers of this small parameter, several nonperturbative effects (in terms of coupling constants) 
may be systematically investigated. 

The ability of the 1/N framework to address the nonperturbative aspects of quantum field 
dynamics has motivated a detailed study of the properties of these systems. In non-equilibrium 
situations, this formalism has been applied to the dynamics of symmetry breaking 11. 38j - 3*51 1401 |4"T] 
and self-consistent semiclassical cosmological models |42( I43| I44[ I45[ I46| I47j . 

In the case of the O(N) invariant theory, in the presence of a nonzero background field (or an 
external gravitational or electromagnetic field interacting with the scalar field) we may distinguish 
the longitudinal quantum fluctuations in the direction of the background field, in field space, 
from the N — 1 transverse (Goldstone or pion) fluctuations perpendicular to it. To first order in 
1/N, the longitudinal fluctuations drop out of the formalism, so we effectively are treating the 
background field as classical. Likewise, quantum fluctuations of the external field are overpowered 
by the fluctuations of the ./V scalar field. In this way, the 1/N framework provides a systematics 
and a quantitative measure of the semiclassical approximation 48 . 

To leading order (LO), the theory reduces to N — 1 linear fields with a time dependent mass, 
which depends on the background field and on the linear fields themselves through a gap equation 
local in time. This depiction of the dynamics agrees both with the Gaussian approximation for the 
density matrix |49[I5()| and with the Hartree approximation [31]. The LO 1/N theory is Hamiltonian 
|34| and time- reversal invariant. It simply does not thermalize. For example, if we set up conditions 
where both the background field and the self-consistent mass are space-time independent, then the 
particle numbers for each fluctuation mode will be conserved. The existence of these conservation 
laws precludes thermalization 51 . 

We note that the failure of the LO approximation to describe thermalization is indicative of a 
more general breakdown of the approximation at later times, where effects of particle interaction 
dominate. Both the distribution of energy among the field modes and the phase relationships (or 
lack thereof) among them affect the way quantum fluctuations react on the background or external 
fields. Therefore, from physical considerations, one can say that a theory which does not describe 
thermalization becomes unreliable for most other purposes as well |52| . 

This is where the next to leading order (NLO) approximation enters. It has been applied to 
quantum mechanics ^ |2] , classical field theory [SHI EU ESI EH| and quantum field theory in 1 space 
dimension [HZj, being contrasted both to exact numerical simulations of these systems, as well as 
against other approximations purporting to go beyond LO. The NLO has been shown to be an 
accurate approximation, even at moderate values of N. 
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The 2PI formalism is also suitable for this question because, provided an auxiliary field is 
cleverly introduced, the 2PI CTP effective action can be found in closed form at each order in 
1/N |1H 158j . We now begin our study with a concrete model. 

3.2 O(N) XX 4 Theory 

To fix the physical ideas, we shall discuss the problem in zero space dimensions, namely, the 
quantum mechanical O(N) model. The system dynamics is described by the Hamiltonian with 
variables Xa and their conjugate momenta P A , where A, B are the O(N) group indices, with 



H 



pBpB + m 2 x b x b + A (x B x B y 



(7) 



The classical action 



S 



1 



A 



dt - { X B X B - M l X B X B - — {X B X B ) 



(8) 



where M5 2 and Aq are the mass parameter and coupling constant. We rescale X B 



Nx B 



N / dt 



x B x B - M 2 x B x B - - (x B x B ) 2 
Discarding a constant term, we may rewrite the classical action as 



(9) 



S = N J dt - I % B ^ B 



M 2 V\ 

7X + ~ XBXB 



(10) 



To set up the 1/N resummation scheme, it is customary to introduce the auxiliary field x, 
writing 




x B x B 



M 2 VA 

7=x + ~ XBXB 



1 2 



M 2 



Vx 



x^^/x 

— + —x B x B 



(11) 



whence 



S = N J dt \ ^x B x B 



X 



M_ x B x B 
X 2 



+ 2A^ 



(12) 



From now on, we consider \ and xb as fundamental variables on equal footing. 

Because of the O(N) symmetry, the symmetric point must be a solution of the equations of 
motion. For simplicity, we shall assume we are within this symmetric phase, and treat X Q ELS cl 
quantum fluctuation. We also split the auxiliary field \ = X + X. i n t° a background field \ and a 
fluctuation field % The action becomes 



S — Sq + Si + S2 + S3 
Sq is just the classical action evaluated at x B = 0, % = X'- 



(13) 



Sn 



N 

~X 



dt 



X 2 ~ M 2 x 



(14) 



Si contains terms linear in x an d can be set to zero by a suitable choice of the background field x : 
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S*2 contains the quadratic terms and yields the tree - level inverse propagators, 



S 2 = N j dt j^i 
Finally S3 contains the bare vertex 

'—TV 



X 

xbxb - -^xbxb + 



— ~ 2 \ 

2X X J 



J d d x {xxbxb} ■ 



(15) 



(16) 



(17) 



To write the 2PI CTP EA we double the degrees of freedom, incorporating a branch label 
a = 1,2 (for simplicity, if not explicitly written, we assume that the label a also contains the 
time branch, i.e., x Aa = x Aa {t a )). We also introduce propagators G Aa,Bb for the path ordered 
expectation values 



G Aa,Bb = ( x Aa x Bb 



(18) 



and F ab for 



F ab = / -a £ 



(19) 



Because of symmetry, it is not necessary to introduce a mixed propagator, for (x a x Bb ) = 0. The 
IP I CTP EA reads 



So 



X 



So 



X' 



^Jdudv ^D Aa ^ Bb (u,v)G Aa ' Bb (u,v) + ^- Cab 5(u,v)F ab (u,t>)} 



-y [Tr InG + Tr In F] + Tq 
where, if the position variable is explicit, c\\ = — C22 = 1, c\2 = C21 = 0, 



(20) 



DAa,Bb(u,v) = N5 A b c ah d 2 x - c abc x c S(u,v), 



(21) 



and c a bc = 1 when all entries are 1, c a b c = —1 when all entries are 2, and c a bc = other- 
wise. When we use the compressed notation, it is understood that c a b = c a (,5(t a ,t(,) and c a fc c = 
c a bc5(ta,tb)5(t a ,t c ). Tq is the sum of all 2PI vacuum bubbles with cubic vertices from S3 and 
propagators G Aa ' Bb and F ab . Observe that Tq is independent of x c - 

Taking variations of the 2PI CTP EA and identifying X 1 = X 2 = X, we find the equations of 
motion 



— i>ABF>ab - -J 



G 



-1 



H II/lo Rh = 

Aa,Bb 2 ' 



N 



ifh 



2X° ab 2 



F 



-1 



ab + 2^ b = ° 



(22) 
(23) 
(24) 
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where D ab (u, v) = c ab [8%. - x («)] S(u, v), 



U Aa,Bb - 2 6 QAa,Bb' U ^ = 2^ (25) 

We shall seek a solution with the structure 

G Aa > Bb = ^6 AB G a \u,v) (26) 
which is consistent with vanishing Noether charges. Then it is convenient to write 

F ab = t_ H ab. U Aa>Bb = 5 AB P ab ; U ab (x,y) = NQ ab (x,y) (27) 
The equations become 



Dab ~ i 
1 

jc ab - i 



G 



-i 



„> + ^* = (28) 



+ Qab = (29) 

ab 



l{x(t)-M 2 }- ±G n (t,t) = (30) 



Observe that 



_2 5Tq n _ 2 5T Q 
ab ~ h 5G ab ' ^ ab ~ h 5H ab 1 j 

These are the exact equations we must solve. The succesive 1/N approximations amount to different 
constitutive relations expressing P ab and Q ab in terms of the propagators. 

The key observation is that in any given Feynman graph each vertex contributes a power of 
N, each internal line a power of N" 1 , and each trace over group indices another power of N. We 
have both G and H internal lines, but the G lines only appear in closed loops. On each loop, the 
number of vertices equals the number of G lines, so there only remains one power of N from the 
single trace over group labels. Therefore the overall power of the graph is the number of G loops 
minus the number of H lines. Now, since we only consider 2PI graphs, there is a minimum number 
of H lines for a given number of G loops. For example, if there are two G loops, they must be 
connected by no less than 3 H lines, and so this graph cannot be higher than NNLO. A graph 
with 3 G loops cannot have less than 5 H lines, and so on. 

We conclude that Tq vanishes at LO, and therefore P ab = Q ab = 0. There is only one NLO 
graph, consisting of a single G loop and a single H line. This graph leads to 

rr = (-^)(-J) ("2^) 22iV (^) 3cabcCde/ / dudvH^MCPiu^Cr'fav) (32) 
Therefore, we get 

Pab = ihc acd c bef H ce G df (33) 

Qab = ^c acd c bef G ce G d f (34) 
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3.3 From correlations to the reduced density matrix 



The 2PIEA yields equations of motion for the (arbitrary time) two-point functions of the theory. 
Given a solution of these equations, in principle we may find the expectation values of a large family 
of composite operators at any given time. Suppose we adopt a coarse-grained description where 
we choose a certain number of these expectation values as the relevant variables to describe the 
system. Then there will be a single density matrix which has maximum von Neumann entropy with 
respect to the class of states reproducing the preferred expectation values. This maximum entropy 
density matrix is the reduced density matrix for the system, and its entropy its correlation entropy. 
The H theorem is the statement that the correlation entropy grows in time, when the correlations 
themselves are evolved using the equations derived from the 2PIEA truncated to some order in 
the 1/N expansion (this statement will be qualified below). 

Let us consider first the simplest case where we describe the system by specifying the values 
of (xbxb) , (xbPb +Pbxb) an d (p B p B ^ (p A = Nx A is the momentum conjugate to xb) at every 
moment of time (by symmetry, we assume (xbxc) = §bc {xbxb) /N, and similarly in the other 
cases). This choice of relevant variables does not utilize or display the full power of the 2PI 
CTP EA, in the sense that the 2PI CTP EA yields equations of motion which, if carried beyond 
LO, account for the build up of non-Gaussian correlations. However, this restriction is assumed 
explicitly or implicitly in most of the work on thermalization, and the focus is placed on the 
shape of the spectrum of the two-point functions at different times. The situation is analogous 
to Boltzmann theory, where the equations of motion can describe the build-up of two particle 
correlation functions from an uncorrelated initial state, but the Boltzmann entropy is defined from 
the one particle distribution function alone. 

At every moment of time we define a maximum entropy density matrix 



p(t) = Z- X e- H * 



where 



and 



H o = 2 {ap B p B + P (xbPb + Pbxb) + lx B x B } 



Z = Tr e~ H ° = z N 
1 



z = 



2 sinh 



\ha 



where 



a 2 = a 7 - (5 2 



(35) 

(36) 

(37) 
(38) 

(39) 



(see Appendix A). The parameter a measures how far the system is from a pure state. To see this, 
observe that 



Trp 2 



Tr e~ 2H ° 
Z 2 



2 sinh 2 



\ha 



sinh [ha] 



N 



tanh 



-ha 



N 



(40) 



So a — > oo yields a pure state, while the state is mixed for any finite a. Therefore, to show the H 
theorem, we must show that da/dt < 0. 



10 



Let us begin by showing that a may be written directly in terms of the expectation values for 
binary products of canonical variables. Indeed, we have 



d 



- 2 — InZ = (p n p 
oa 



BR 



Nh 7 



tanh 



2a 



v Nh a 

- 2— InZ = {x b xb) = 

°1 tanh 



\ha 



d 



- 2— In Z = (x B PB + Pbxb) = ~ 
OP tanh 



2a 

Nh (3 



\ha 



Therefore 



a 



(41) 
(42) 
(43) 



4 (p B p B j {x b xb) - (xbPb + Pbxb) 
As an aside, these formulae show that 

(Ho) 

therefore the von Neumann entropy is 



Nh 



tanh 



:ha 



Nha 



2 tanh 



\fia 



S = - (hip) = (flo) + lnZ 



Nha 



2 tanh 



\ha 



AT In 



sinh 



-ha 
2 



iVln2 



(44) 



(45) 



(46) 



Observe that 



dS 



Nh 2 a 



\ha 



< 



(47) 



da 4 sinh 2 

So again, to obtain a H theorem we must show that a is non-increasing in time. 
3.4 The physical basis of the H theorem 

It is clear that if we could solve the exact evolution, the H theorem would be manifest: Given p (t ) 
at the initial time to, solve the exact Liouville equation up to a time t\. Let p(t\) be the result. 
We extract the new expectation values from p{t\) and use them to construct the new maximum 
entropy density matrix p (ti). Then S [p (t\)] > S [p (ti)] , by definition, and S [p (ti)] = S [p (to)], 
because the exact evolution is unitary, thus the H theorem. 

What we need to determine is whether, given the approximate dynamics for the expectation 
value provided by the 1/JV scheme, the H theorem still obtain. In fact, we know to LO it does 
not. From the angle of addressing the thermalization issue, the LO approximation totally misses 
the point. 

Physically, we expect to obtain a H theorem because our description of the system is incomplete, 
in that it ignores higher correlations. We may rightfully call this entropy correlation entropy, as 
explained in the Introduction. The system as described by a finite order 1/N approximation is an 
open system, in the following sense: Suppose we try to reproduce it in the lab. In order to have 
the different correlations evolving according to the 2PI 1/N equations appropriate to the desired 
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order, rather than the full Schwinger - Dyson hierarchy, we must keep nudging it to conform to 
this artificially created condition (again, due to our inability to comprehend the complete picture). 
So there is energy and (in principle) entropy flow in and out of the system besides correlation 
entropy production. The sum of the entropy produced in the system, and that exchanged with the 
environment (which keeps the system on course) need not be positive. 

In the LO case, as the approximate evolution (described by a quantum Vlasov equation [55] ) 
is unitary, there is no net entropy production and no //-theorem. This implies that the two 
sources of entropy change must cancel exactly. Note that in such a case, if one works with a Fock 
representation, for boson fields, the number of particles increases with time and can be used as 
a measure of field entropy (what Hu and Pavon [3U| called an intrinsic measure of field entropy, 
as described in Sec. IIA above. See also (HI]-) However, this should not be confused with the 
correlation entropy of the Boltzmann kind under study for which the //-theorem is defined. 

Of course, we expect better approximations to be closer to the exact dynamics, therefore re- 
quiring less external control of (or rather, tampering with) the system. This reduces the entropy 
loss to the environment. Eventually correlation entropy production becomes the dominant factor, 
and a H theorem is obtained. 

To summarize, the relevant question is not whether there is a H theorem for a given choice of 
relevant variables, which is obvious, but rather if the NLO approximation is good enough to make 
it manifest, or we must go even higher. 



3.5 Statement of the H theorem 

To continue our investigation, first note that 



A = {4 (p B p B ^ {x B x B ) - {x B PB + Pb^b) 2 ] 



2*2 



da 



tanh 



ha I sinh 



\ha 



dt 



(48) 



So we must show that the left hand side is positive. Since the 1/N scheme proceeds via taking 
the expectation values of the canonical equations of motion, we may use them to simplify this 
expression. The canonical equations are (for simplicity, we use the form without the auxiliary field) 



Therefore 



P A 



-N 



■A P A 

x =lv 



A 



M + - (x B x B ) 



(49) 
(50) 



d_ 

dt 



dt 



(x B x B ) 



N 



{XBPB + PBXb) 



p B p B 



-N 



M 2 (x B PB +Pbxb) 



A 



Recall the canonical commutation relations 

1 



xbPb = 2 {xbPb+Pbxb) + "2 



{(x C X C ) XBPB + PBXB (x C X C )) 



ifi 



(51) 
(52) 

(53) 
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so 



Finally 



so 



1 ifi 
PBXB = 2 (xbPb+PbXb) ~ y, 



-iV M 2 (xbpb + Pbxb) 



+- {{xcxc) (zbPb +Pbxb) + (xbPb +Pbxb) (xcxc)) 



j t {xbPb+Pbxb) = (pV) - 2N 



M 2 {xbxb) + 77 ((xcxc) 



A = — {4^p B p B ^ (x B x B ) - (xbPb +Pbxb} 2 } 

= XN ^2(x B pB +Pbxb) ({xcxc) 2 ) 

- (xbxb) ((xcxc) (xbPb +Pbxb) + (xbPb +Pbxb) (xcxc))} 

This may be rewritten in the following suggestive way 



A |2 ((x c x c ) 2 ) ^ (x B x B ) - (xbxb) ^ ((xcxc) 2 )~j 
(xcxcf 



(-A) (x B x B 



so the H theorem reads 



Or, equivalently, 



dt 



(x B x B y 



d_ 
~dt 



xcxc) 



(x B x B y 



< o. 



dt 



(x c x c ) 2 ) - (x B x B ) 2 



< 



(x B x B ) 

Observe that this expression is invariant under a rescaling of the field. 



4 Correlation entropy production in weakly coupled theories 

The results from the last section may be summarized as (cfr. Eqs. ()47|) . (|48|) and 1)591)) 



dS f-X\ 

n = {w ](JtaQh 



(x B x B ) - 



(xcxc)' 

(x B x B y 
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We shall not attempt to give a general proof that this quantity is non-negative, but only check that 
it is so in some simple cases. 

We must first express the expectation value in Eq. (|62l) to the required order in 1/JV. Recall 
that the q— number auxiliary field x was introduced as a formal Gaussian process with correlation 
functions 

( x (t)) x = M 2 + ^x B x B (t) (63) 

(x (t) X (t')) x = (M 2 + \x B x B (t)) (M 2 + \x B x B (0) +™6{t- t>) (64) 

where () denotes an expectation value with respect to the x variable alone. Taking a further 
expectation value over the quantum state of the x fields, we recover Eq. (|3fl|) and 

^H u (t, t') = ^ [{x B x B (t) x B x B (t')) - (x B x B (t)) (x B x B (*')>] + (t - t') (65) 
Therefore the H theorem boils down to showing that 

«l[«"M)-iAi(0)] < „ 

Mi (x B x B (t)f - ' 

It is clear that there is no correlation entropy production in free theories, so if the H theorem 
is to hold as a categorical relationship, then it must hold already to lowest order in the coupling 
constant. Therefore we may expand H 11 (t, t) in powers of A, keeping only the lowest nontrivial 
order. Iterating the equation of motion 

H ac (t, t') = i\c ac 5 (t-t')-xj ds c ab Q bd (t, s) H dc (s, t') (67) 



we get 



where 



H u (t, t) = iX5 (0) - i\ 2 Q u {t, t) + i^-y (68) 



dt' 



(Qu (t,t')) 2 - (Qu (t,t')) 2 

4 / „ , „\4 



(G 12 (t,t')) -(G 11 (t,t) 



(69) 



We shall investigate the H theorem in two limiting cases, first at early times for an arbitrary 
(diagonal in occupation number) initial state, and then at late time for vacuum initial conditions. 

4.1 The H theorem at early times 

To investigate the meaning of this expression, let us expand the LO propagators in terms of mode 
functions 

/(t) = J=e--'. (70) 
V 2lo 

Concretely, if the initial state is diagonal in occupation number, we have 
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G 11 (t, t') = J_ e -^(*>-*<) + - C o S co (t - t') 



n 



G 12 (t, t') = J- e ^(*-f ) + - cos a; (t - t') , 
2uj uj 



where n = (0)0). Therefore 



(See Appendix B). Since 



2n 3 + 3n 2 - n 



sin 2a;t + 



(x B x B ) = 



2 (2n 3 + 3n 2 ) + 4n + 1 



sin 4wi 



1 + 2n 
2w 



(71) 
(72) 

(73) 
(74) 

(75) 



is constant, the relevant derivative is (i)d^f/dt, which is proportional to d(— J tpdt) jdt. The H 
theorem [(i) d^/dt < 0] will obtain if (d/dt) f ipdt > 0. Now, since ip (0) = 



(76) 



which is indeed positive at short times. 



4.2 The H theorem at late times 

Let us assume vacuum initial conditions. The lowest order propagators may be written in terms of 
mode functions 



G (t,t')=F(t > )F* (t<) 

G 12 (t, f) = F* (t) F (t') . 
* = J* dt' [f a (t f ) F* A (t) - F A (t) F* A (t') } . 



We have 



{x B x B ) = h\F(t)f 



So the relevant inequality is 



.d_ 1_ 

''di\F(t)\' 



■ J* dt' [f a {t') F* A (t) - F A (t) F* A (t') } < 



Or, equivalently, 



Now 



dt Jo 1 F A (t) F* 4 (t) J " 



(77) 

(78) 
(79) 

(80) 

(81) 

(82) 
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J = - Al Io dt ' { F ' {t) l^§- F *' {t) 



F* 4 (t') 
F* 5 (t) 

The mode functions may be expanded in terms of WKB modes 



(83) 



F(t) = a(t)f(t) + f3(t)f*(t) (84) 

F'(t) = a(t)f'(t)+P(t)r'(t) (85) 
At long times, the integrals will be dominated by the non-oscillatory terms 

F*(t')~6a 2 (t)f3 2 (t)\f(t)\ 4 (86) 

J ~ \f (*)! 4 \ F ' (*) F * 5 (*) « 2 (<) ^ (*) " F *' (*) (*) «* 2 (*) (*)} • (87) 

F Wl 

We make the analogous approximation 

F' (t) F* 5 (t) = a (t) f (t) F* 5 (t) + (t) /*' (t) F* 5 (t) 

~ 10 |/ (t)| 4 a* 2 (t) (3* 2 (t) {\a (t)\ 2 f (t) f* (t) + |/3 (t)\ 2 /*' (t) / (t)} (88) 



and get 

J <■ 

\I 

But 



J ~ \f (t)\ 8 \a (t)| 4 |/3 (i)| 4 {\a (t)\ 2 - |/3 (i)| 2 } {/' (t) /* (t) - /*' (t) / (t)} . (89) 



|a(t)| 2 -|/3(t)| 2 = l (90) 

f(t)r(t)-r'(t)f(t) = -i. (9i) 



So 



J ~TT^ol/(*)| 8 |«(*)| 4 |^(*)| 4 <0 (92) 
F Wl 

The existence of a H-theorem at the NLO is as reassuring as the absence thereof at the LO. 
It is also interesting to observe that entropy production is (formally) a quantity of the order of 
iV -2 ; a strict (formal) 1/N expansion would neglect such a quantity, but it would be generated in a 
numerical solution of the full 1/N equations of motion. Therefore, it is suggestive that the statisti- 
cal mechanical description of an interacting quantum field at the NLO level of approximation to a 
2PI CTP effective action may permit "thermalization" , as claimed by |211 122j based on numerical 
results. Note that our analysis corroborating this claim is predicated upon the restrictive conditions 
discussed above. Extension of this work to full field theory in conjuction with the reported results 
of ^3] is under investigation. 
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5 Appendix A: Computing the partition function 

To compute Z, it is best to diagonalize H . Introduce new canonical variables 

£a = (cosh u) pa + (sinh u) xa 



(93) 



so that [£a,Vb] = \pa,x b ] 



rj a = (sinh it) pa + (cosh u) xa 
(—ih)5AB- Then 

PA = (cosh u) ^a — (sinh u) tja 



(94) 



(95) 



and 



where 



xa = — (sinh u) £a + (cosh u) tja (96) 



H = (~) {A£ b £ b + B ( VB ^ B + £ B 7?b) + Ct,bVb} (97) 



A = a cosh 2 u + 7 sinh 2 n — 2/3 cosh u sinh u (98) 

B = — (7 + a) cosh u sinh u + (3 (cosh 2 u + sinh 2 uj (99) 

C = a sinh 2 u + 7 cosh 2 u — 2/3 cosh u sinh u (100) 
To diagonalize the Hamiltonian, we require B = 0, namely 

20 

tanh2u = — — (101) 
7 + a 

Now Z is the usual partition function for N harmonic oscillators with frequency uj 2 = C/A at 
inverse temperature A. 



6 Appendix B: ijj (t) 



*(*) 



dt' 



e iu(t-tf) + 2n cos u(t-t') 



-Mt-t') + 2n cos uj(t-t') 



16w 4 Jo 

/ dt' (32n 3 cos 3 w (t - t') sin w (t - t') + 24n 2 cos 2 u (t - t') sin 2uj (t - t') 
8a; 4 jo 

+8n cos uj {t — t') sin 3a; (t — t') + sin 4a; (t — t') } 
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ft 

J dt' 1 8 (2n 3 + 3n 2 ) cos 2 w (i - t') sin 2a; (t - t') 
+8n cos oj (t — t) sin 3a; (t — t) + sin 4a; (i — i') } 



i 

8^4 



y tit' 1 2 (2n 3 + 3n 2 ) [2 sin 2a; (t - t') + sin 4a; (i - i') ] 



+4n [sin 4a; (t - t') - sin 2a; (t - t')] + sin 4a; (t - t')} 



i 

8w4 



dt' 1 4 sin 2a; (t - t') 



2n 3 + 3n 2 



n 



+ sin 4a; 



(t - t') [2 ( 



2n 3 + 3n 2 ) + An -f{l|i) 
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